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We present designs for multipole ion traps based on a set of planar, annular, concen- 
tric electrodes which require only rf potentials to confine ions. We illustrate the desirable 
properties of the traps by considering a few simple cases of confined ions. We predict that 
mm-scale surface traps may have trap depths as high as tens of electron volts, or micromo- 
tion amplitudes in a 2-D ion crystal as low as tens of nanometers, when parameters of a 
magnitude common in the field are chosen. Several example traps are studied, and the scal- 
ing of those properties with voltage, frequency, and trap scale, for small numbers of ions, 
is derived. In addition, ions with very high charge-to-mass ratios may be confined in the 
trap, and species of very different charge-to-mass ratios may be simultaneously confined. 
Applications of these traps include quantum information science, frequency metrology, and 
cold ion-atom collisions. 



I. INTRODUCTION 



Ion traps are ubiquitous tools in physics and chemistry, with applications including quantum 
information science [1J, precision measurement [|2l|3l, and mass spectrometry Q. The two main 
classes of ion trap are Paul traps, which use an oscillating electric field, usually combined with a 
static electric field, and Penning traps, which combine a static magnetic field with a static electric 
field. In many applications, such as quantum information science, the ions are laser-cooled, which 
can cause them to condense into an ordered crystal. Even at low temperature, however, the ions in 
both types of trap suffer an unwanted excess motion: rf-driven micromotion, in the case of Paul 
traps, and rotation of the crystal due to the crossed electric and magnetic fields in Penning traps. 
In linear Paul traps, the ions form a one-dimensional chain, and micromotion along the axis of 
the chain is minimal, but this is not so for crystals of two and higher dimensions. So far, there 
has been no method of preparing a 2D or 3D crystal of motionless charged particles in a single 
trapping region. 



One approach to building a trap in which ions are largely free from excess motion is a so-called 
multipole ion trap. The most common variety of such, the linear multipole trap, is made of several 
parallel rods, the voltage on each being either 180° out of phase with its neighbor or grounded. 
Although laser-cooled crystals have been observed in such linear multipole traps [5J, one encoun- 
ters the unwanted effect of the total potential minimum being displaced from the minimum of the 
rf pseudopotential by the static axial confining field [6j. This leads to greater micromotion than 
would otherwise exist. Two particular deleterious effects of this micromotion are errors related to 
the Doppler shift in atomic clocks based on trapped ions [6], and a fundamental upper limit to the 
temperatures of ion-atom mixtures Q. A recent report shows that micromotion- free parallel ion 
strings can be formed with additional rf potentials 10. However, we are interested in the question 
of whether it is possible to create a multipole trap potential with a minimum on the axis of sym- 
metry, and in which all ions occupy a single trapping volume. We refer to this situation as an ideal 
multipole ion trap. 

Our solution to this problem is inspired by novel ion trap designs based on electrodes that lie 
in a single plane |[9l-PTTT|. These surface-electrode or planar traps can be microfabricated using 
standard optical lithography. The advantages of miniaturization include higher interaction rates, 
integration of control fields into the trap structure [172141411 . and the possibility of moving ions 
between memory and processing zones [15]. A layered planar multipole trap has been developed 
lfT6l ; however, its linear geometry implies that axial confinement is still necessary, both leading 
to rf-driven motion and impeding optical access. We therefore ask the question of whether a 
surface-electrode multipole trap (SEMT) can be built. 

The remainder of this article is organized as follows. In Sec. |n} we outline the general mathe- 
matical approach to the problem. In Sec.linJ we present solutions obtained in a number of ways 



for SEMT's. Sec. IV is devoted to understanding the properties of the traps, in particular relating 
experimentally adjustable parameters such as the drive frequencies and voltages to properties such 
as the trap depth and micromotion amplitude for a particular ion. We conclude with Sec.|V} which 
contains a discussion of the possible applications of these traps. 
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FIG. 1: Our model trap consists of a set of planar, annular electrodes. The cylindrical coordinate r is 
measured from the center of this structure, and the coordinate z measures distance above this structure. The 
electrode widths are arbitrary and are specified here using r\, r%, etc., as shown, which is possible in the 
gapless approximation which we use. The voltages on each electrode V\, Vi, ... may be set to any value; 
the "ellipses" indicate that an arbitrary number of rings may be used. The ring electrodes are assumed to 
be surrounded by an infinite grounded (V = 0) plane. 



This work deals exclusively with traps based on concentric rings. This choice is based on the 
fact that a ring geometry allows one to trap ions in three dimensions using only rf (and no nonzero 



either an rf voltage or zero voltage is applied. We assume that the distance between electrodes is 
zero and that the outer electrode is surrounded by a plane of infinite extent. Much work has been 
published on analytical solutions for the electric potential of surface-electrode ion traps IfTTT - ETII : 
here, we find it convenient to use the form of the solutions presented in Ref. IfTTl . 

We first review the general solution for our trap geometry. Employing cylindrical coordinates 
z and r, where the origin for z is on the trap surface and for r is on the axis of symmetry, the 
electric potential is written 



where Jo is the Bessel function of zeroth order and A (k) = J2i=i The A, are given by 



II. GENERAL APPROACH 



dc) potentials IfTTl . Our model trap (Fig. [TJ consists of a set of fiat conducting rings to which 




(1) 



Ai(k) = Vl (r i+ iJi(kr i+1 ) - nJ^kri)) , 



(2) 
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where Vi is the amplitude of the rf voltage applied to the i th electrode, J\ is the Bessel function of 
first order, and the are the radii describing the trap geometry (as shown in Fig. [I]). In our calcu- 
lations, we make the pseudopotential approximation, in which the time-averaged pseudopotential 
is given by 



Q_ 2 



^{^r) = -^—\E{z,r)\\ (3) 



where Q is the ion's charge, m is the ion's mass, Vt is the frequency of the rf potentials applied to 
the electrodes, and E is the electric field produced by the trap electrodes. This approximation is 
valid when the time scale for the variation of the ion's position and velocity is much greater than 
the rf period [J22]]. 

We define the variable y as y = z — z , where z is the chosen height of the trap center above 
the electrode plane. The trap center then is located at the point (r = 0,y = 0). The potential may 
be expanded in a power series about the trap center as 



$(z, r) = (cyo r o + c ylr0 y + c y2r0 y 2 H ) + r 2 (c y0r2 + c ylr2 y + c y2r2 y 2 H ) 

+ r 4 (c y0 r4 + c y i rA y + c y2 r4y 2 H ) H • (4) 

Our goal is to construct a pseudopotential, the lowest order of which in r and y scales as r p and 
y** for some p, ( > 2. By Eq. p and ( must be even. This is tantamount to setting the terms 
c y 2ro an d c y0r2 to zero. We do not venture to also null the coefficients on the "cross terms" which 
involve products of powers of r and y. 

Understanding the relationship between r and y is important because $ has no closed-form 
solution, except when r = ifTTll . Let us see how we can use derivatives along y to null the 
potential to some order along r. We take successively the gradient and then the divergence of 
Eq. evaluating it at the trap center. Noting that the Laplace condition states V 2 $ = for all r 
and y, and that any term containing r or y to any power goes to zero at the trap center, one arrives 
at the following relationships: 
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— 4q,0r2 + 2c y 2r0; 

= 64c y 0r4 + 24cy 4r0 ... 



(5) 



and so forth. We note that the form of the divergence in cylindrical coordinates implies that 
CyOri < Cyiro for any i. These conditions assure that the cross-sections along r and y are super- 
quadratic, and that the lowest-order nonvanishing term in the expansion of $ is at least of cubic 
order (in r or y individually, or in a product of coordinates). We aim therefore to find solutions to 
the set of equations of the form 



<9 fc $ 



z=z 



0, (6) 



dz h 

where k ranges from 1 to the desired order n. The k — 1 condition implies that the electric field 
is zero at the trap center. As an example, if we wish to form the lowest-order multipole trap, we 
must null the k — 1 and k = 2 derivatives. To null the pseudopotential along both directions to a 
given order, all derivatives must go to zero. However, to null only along r, for instance if one is 
satisfied with a quartic trap along y, only the even derivatives greater than k — 1 must be zero. 



III. CONSTRUCTING MULTIPOLE TRAPS 

The equations of Eq. ([6]) are satisfied by a correct combination of the electrode widths r, ; or the 
voltages Vi. According to Eq. ([T]), the potential is linear in each applied rf voltage; however, it is 
quite nonlinear in y (or z) and in the r^. Therefore, we first solve the system for the Vi after fixing 
the electrode widths. We shall then examine the alternative approach. 



A. Varying the voltages 

As a first example, let us construct a lowest-order multipole trap, meaning that the quadratic 
terms in r and y, but not necessarily the cross-terms depending on products of powers of r and 
y, go to zero. In addition to setting the first and second derivatives of \l/ with respect to y to zero 
at y — (z — z ), we are free to choose the value of the potential at some point; therefore we 
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TABLE I: Table of solutions to the system of Eq. ([6]), with the potential nulled along both y and r. k m i n 
represents the lowest order term in the expansion of &(y) that is nonzero. The lowest-order nonvanishing 
terms in r and y are given by a and (, respectively. The number of voltages needed to satisfy the increasing 
number of conditions increases linearly with k m i n . These solutions are only valid for the case in which 
each electrode has unit width (and the center electrode has unit radius), V\ = 1, and zo = 1. 

specify V\, the voltage of the center circular electrode: $(2; = 0) = V\, In all, we seek to satisfy 
the following conditions: 



$1 

\y= 


-Z0 






dy 


y=( 


<9 2 $ 




dy 2 


y=o 



= v 1 

= 



We first choose a simple layout in which each electrode has the same width, equal to the radius 
of the center electrode, chosen to be 1. We also choose V% = 1 and z = 1. Satisfying Eq. ([7]) 
requires solving for two unknowns: the voltages V 2 and V 3 . Taking the derivatives of $ and 
forming the system, we arrive at the solution V 2 = —0.29, V 3 = 4.82. We observe that indeed the 
lowest-order nonvanishing term in both and ^(y) is proportional to y 4 . 

Our next step is to calculate the pseudopotential along r, which cannot be solved analytically. 
However, for small r, \I/(r) is well approximated by a polynomial containing one or two terms. 
Therefore, we calculate a set of points using numerical integration, and fit to a function containing 
those terms. In Table I, we provide the solution for traps of increasing order in k, and in Fig. [2j 
we plot ^(z) and \l/(s) for the lowest-order multipole trap. 

The voltages in Table [j] increase quite rapidly with increasing order k. Moreover, one may 
not be interested in nulling the potential along z, but only along r, for instance if one wishes to 
study planar ion crystals for which the confinement along z is much tighter than along r. We 
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FIG. 2: Plots of ^(z) at r = (top) and \t(r) at z = zq (bottom). These plots are scaled using experimental 
parameters given in Sec. IV In particular, V\ = 100 V and the drive frequency Q,/(2ir) = 2.86 MHz. 
Although ^(r) is well approximated for r <C 0.001 m by a function proportional to r 6 , this fit which 
contains the edges of the trap contains terms up to r 14 . 
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TABLE II: Table of solutions to the system of Eq. Q, with the potential deliberately nulled only along r. 
Here, k m - m represents the lowest order nonvanishing even term in the expansion of ^(y). The lowest-order 
term altogether has k = 3. The other symbols and parameters are as defined in Table [I] 



now present solutions for the correct voltages with only even orders nulled. Since the coefficient 
c y 3ro is not zero, we expect ^ to have a leading order in y of £ = 4, even though by nulling the 
even derivatives, we can increase the leading order p in r to an arbitrarily high number. Table [TT] 
provides a list of these solutions. We refer to the lowest-order trap described in Tables I and II as 
Trap A. 

There are numerous combinations of electrode widths, z values, and electrode voltages that 
will result in a multipole trap of lowest order. For example, it is possible, at least for Trap A, to 
vary z so that one of the applied voltages approaches zero. From an experimental point of view, 
it may be greatly preferred to limit the number of different rf voltages applied. For the same trap 
geometry as Trap A, but with an ion height of 0.88, the solution is V% = 1.2 x 10~ 4 , V 3 = 4.76. 
Also, there is no reason that all electrodes must be of the same width. Indeed, in the next section 
we shall see that two electrodes of unequal width can carry the same rf voltage, also resulting in 
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a multipole trap. Let us consider, as another example, a situation in which r% — 1, r 2 = 3, and 
r 3 = 5. For z = 1.7265 and V x = 1, a value of V 3 = 2.61787 results in V 2 = 0. 

B. Varying electrode widths 

Although the most direct and convenient route to a solution is to fix the electrode widths and 
permit the voltages to vary, from an experimental point of view it may be preferred to build the 
correct field curvature into the electrode structure, so that a single rf voltage (or ground) may 
be applied to all electrodes. As noted above, this problem is more difficult than solving a linear 
system for the unknown voltages, but we demonstrate here that it is possible, at least for a small 
number of electrodes. 

Let us construct a trap with $(y) oc y 3 in the lowest order. We assign, again, a unit width to the 
inner electrode (r x = 1). We set the potential of the first concentric ring to V 2 = 0, guessing that 
this is possible based on the result for Trap B above. The potential of the outer ring is V3 = V\ — 1. 
Since we are still free to choose z , there are an infinite number of pairs (r 2 ,r 3 ) that will work. 
One possible solution is (r 2 = 3.49, r 3 = 8.63), which gives zq = 2.46. Again, we observe from 
fits to the pseudopotential that p = 6 and £ = 4. Henceforth, we refer to this trap as Trap B. 
Although this approach does appear to work, we do not venture to null higher orders. 

Although, from an experimental point of view, it is inconvenient to use multiple rf signals 
to drive a trap, the trap fabrication will never be perfect, and there will be deviations from the 
solution found here due to the gapless and infinite-plane approximations [20]. Therefore, we wish 
to demonstrate that adjustment of the rf voltages can compensate for imperfections in the electrode 
geometry. We do this using a trial-and-error approach simply as a proof of concept. Let us set r 3 
to 8.8, rather than 8.63, and vary V 3 . We observe that V 3 ~ 0.97 leads again to the desired trap 
curvature. 



C. A higher-order trap 



We provide one further example. Higher-order traps were derived in Sec. Ill A, but here we 
again "guess" that the electrode widths should continue to increase further from the center in order 
to limit the increase in the applied voltages with increasing r. We then solve for the voltages. The 
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TABLE III: This table summarizes the five example traps that were obtained in Sec. Ill All traps are 
lowest-order multipole traps except for Trap C. 



set of radii used is 1, 3, 6, 10. Again desiring one voltage to go to zero, we set the ion height to be 
2.69625, yielding a solution V 2 = 1.28, V 3 « 0, V A = 4.50. We refer to this as Trap C. As a guide 



for the rest of the paper, we provide in Table III a summary of the three example traps. 

We conclude this section by noting that the optimization of electrode shapes for construction 
of a multipole trap may be done in a more sophisticated way using the approach of Schmied et al. 
Il23l . For the lowest-order trap, the results of such calculations agree with ours E6l . Finally, since 
we have seen that imperfections in the trap dimensions can be compensated by adjustments in 
voltage, it may be possible to devise a scheme which varies electrode widths and voltages, perhaps 
optimizing for trap depth, as well as curvature. The complication is that one may not necessarily 
want simply to optimize for trap depth; nulling of micromotion may be more important. The latter 
depends not only on the electrode widths and voltages, but also on the number and types of ions, 
and the drive frequency. The next section is devoted to these considerations. In particular, we 



perform an approximate optimization of the lowest-order trap in Sec. IV B 



IV. PROPERTIES OF THE TRAPS 



So far, we have shown a number of ways in which surface-electrode multipole traps can be 
constructed. We now consider the properties of these traps given realistic experimental parame- 
ters. Of interest are the relationships between the overall scale of the V$, the rf drive frequency VL, 
the trap depth D, and the micromotion amplitude of a particular ion, A. Because these traps are 
highly (and deliberately) anharmonic, we do not consider the "secular frequencies" of the ions, 
as these depend, in a multipole trap, on the kinetic energy (or temperature) of the ions and on the 
number of ions trapped. The angular modes of a circular Coulomb cluster, which can be formed 
in SEMT's, are described in Ref. Il24ll . 

In a quadrupole trap, the stability parameters (commonly written "q" and "a") are constant as 
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long as the ion is within the quadratic region of the trap. Since the subject of this article is ion 
traps that require no dc voltages for confinement, we focus here on the parameter q. Although 
there are numerous regions in which ions have stable trajectories, the trap is typically operated 
with a q < 0.3. By contrast, with all multipole traps, one cannot simply solve for a value of f2 
based on fixed formulas for trap stability. Instead, for each value of f2, one must first compute the 
positions of the ions in the trap. The reason is that the parameter that quantifies the stability of the 
ions' trajectories depends on the positions of the ions. For general multipole traps, a quantity r] 
which quantifies local adiabaticity is defined as 



2Q\VE\ 



V = ^7^, (8) 



where E is the magnitude of the (local) electric field [22] . The values of 77 generally considered to 
be stable lie between and 0.34, in close analogy to the parameter q. The micromotion amplitude 
A for a particular ion is given by 



A 



QE (r ) 



(9) 



mf2 2 

where E (r ) is the magnitude of the electric field oriented along f at the ion location [22 J. We 
approximate that the magnitude of the field along z is negligible. This is because, for small 
crystals, the ions will lie in a plane parallel to the trap electrodes, due to the stronger confinement 
along z. 

A very simple, analytical analysis of ion behavior close to the trap center is possible with 
SEMT's. The reason is that, with only rf voltages required for stable ion confinement, the pseu- 
dopotential at (y = 0) expanded about r = may be approximated by a single term. In the case 
of the lowest-order traps, that term is proportional to r 6 . We wish to incorporate the length scale 
of the trap, represented by the ion height z , into the calculation. The potential $ also depends 
linearly on the voltage scale, which we allow to be set by V\. Let us then use this expression: 

$( r ) = ^l r 4 . (10) 

Here, c r is a dimensionless constant that contains information about the effects of all the electrodes 
on the trap curvature, as well as the ratios of their voltages. Depending on the parameter that 
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defines the voltage scale (for example V\, or the maximal voltage on any electrode) and on how 
one wishes to define the length scale (for instance, using r\ rather than z ), the meaning and value 
of c r will change. 

It is now possible to find simple expressions for the other pertinent quantities, including the 
electric field, its derivative, and the pseudopotential. The pseudopotential ^(r) at y = is given 
by 

Ac 2 2 V 2 

and the pseudo-force T(r) = —d^/dr given by 

> 24c 2 r Q 2 V 2 5 

T r = — m~r- r • (12) 

We will use the quantity T when discussing multiple confined ions. The parameter i] at some 
value of r may be easily computed from Eq. ([8]), while the electric field may be gotten from 
differentiating Eq. ( flO| ), and A obtained from that result and from Eq. Q. 

This basic theory can be applied readily to two simple, but important cases: a single trapped ion 
with some nonzero total mechanical energy, and two ions at absolute zero. In the remainder of this 
section, we consider these two cases, as well as crystals consisting of small numbers of ions. We 
assume the crystallized ions have well-defined positions, despite the cylindrical symmetry of the 
trap, as was observed in [17J. These simple analyses shed light on some of the special properties 
of multipole traps, and especially SEMT's, which lend themselves easily to a simple description. 



A. A single ion 



The region of stability for a single ion in a quadrupole trap is strongly limited by the stability 
parameter q. The basic scaling laws are as follows. While q oc QV/Amr 2 Vt 2 , the trap depth D 
goes as D oc qV. Here, V is the single applied rf voltage and r is a constant with dimensions 
of length that sets the scale of the trap; it may also correspond to a specific dimension (such as 
the distance of the ion from a particular trap electrode). Often, one sets V to the maximum value 
attainable in one's experiment (to maximize trap depth) and then selects an so that the resulting 
trap is stable. This effectively sets an upper limit on the ratio Q/m that can be confined in the 
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trap. 

The situation in an ideal multipole trap is very different. The parameter 77 increases with 
distance from the center of the trap and is actually zero at the trap center. For a single ion, then, 
the upper limit on Q/m is determined by its distance from the trap center, due to its total energy 
or (at absolute zero) stray fields. We can derive a set of scaling laws that describe the properties 
of a trapped ion. Suppose that the ion has a total mechanical energy £. Then its maximum extent 
can be obtained by simply setting E = ^(r max ): 



max \4c 2 Q 2 V 2 
Invoking Eq. ([8]), together with Eq. ( 10), we obtain an expression for 77: 



1/3 



-1 = 24(^1 • 04) 



As an example, let us consider some concrete numbers. The value of c r for Trap A is 0.038. Let 
us assume values of Zq — 1 mm, f2 = 2ir x 10 6 s _1 , and £ = k B x 1CT 3 J (for an ion laser-cooled 
to around 1 mK). For a mass of 100 amu, the charge that could be held on the ion, for 77 = 0.3, is 
Q = 4.0 x 10 4 e c , where e c is the elementary charge. It should be pointed out that an ion with any 
nonzero value of Q less than this will also be confined, although the trap depth will fall as Q 2 . It 
is not difficult to calculate 77 for a single ion with any trap order n. It is 

re — 2 

_ 2Qn 2 cV UmQ 2 z 2n S\ 2 - 2 
V ~ mWz™ V Q 2 n 2 c 2 V 2 ) ' 

It is also possible to compute a value for A for the single ion with nonzero energy. Using 

Eq. with a little algebra, we find 

2 IT 

A -Il1Vn <16) 
This expression does not depend on the charge Q. Therefore, there is no disadvantage from the 
standpoint of micromotion to storing a highly-charged ion. This expression also holds for all 



orders of traps, while e.g. Eq. (14) holds only for the lowest-order traps. Remarkably, A also 
does not depend on details of the electric potential such as c r and V\. The implication is that the 
lowest-order multipole trap is as good as a higher-order trap, from the standpoint of micromotion 
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FIG. 3: Logarithmic plots of the trap depth D (left) and micromotion amplitude A (right), as functions of 
for two 88 Sr + ions confined in Trap A with zq = r\ = 1 mm and V\ = 100 V. The trap depth scales as 
D oc Q~ 2 , and the amplitudes as A oc S1~ L15 . 

reduction, but only for a single ion. 

B. Two cold ions 

Another case that lends itself easily to analysis is the confinement of two ions at absolute zero. 
The ions, by Coulomb repulsion, are each located some distance r eq from r = 0. We assume, 
again, that the two ions lie in the plane parallel to the trap electrodes. This is approximately 
true in a real trap. We begin with some examples designed to elucidate the unique properties of 
SEMT's, and then present the general theory. 
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We first compute values of Q, A, and D for Trap A. Rather than an exhaustive search over 
all possible parameters, we examine a few special cases. In particular, we consider two length 
scales: r\ — 1 mm and r 1 = 0.1 mm. We choose a voltage scale V\ = 100 V, resulting in 
V2 = —29 V and V3 = 482 V. Finally, we assign to the ion the charge and mass values of 88 Sr + . 
Plots of D and A as functions of are given in Fig. [3j Unsurprisingly, the trap depths follow 
a simple proportionality to 1/fi 2 , according to Eq. The amplitudes follow a different power 
law, roughly 15 , which is explained below. It is most important to note that the highest depth 
corresponds to the highest level of micromotion, and vice-versa. We can also identify a reason 
why one might wish to scale these traps down. With z = 1 mm, an appropriate drive frequency 
is 1.5 MHz, leading to A = 148 nm. With zq = 0.1 mm, by contrast, a drive frequency of 15 MHz 
(which leads to the same trap depth) gives A = 39 nm. It is useful to compare this to the situation 
in a quadrupole trap, in which A « qr eq /2, where r eq is the distance from the rf null, and q ps 0.3 
is a typical stability parameter. Let us introduce a relative micromotion amplitude K = A/r eq . The 
idea is to quantify the level of micromotion, taking into account the overall size of the ion crystal. 
For quadrupole traps, then, K ps 0.15 (in the radial direction, in the case of a linear trap) under 
typical conditions. For two ions in a quadrupole trap, a typical value of the distance of each from 
the trap center is r = 2 /mi; in this case, A = 300 nm. Let us choose a small, but experimentally 
feasible value for the trap depth: D = 100 meV. The frequency that leads to this depth, for Trap A, 
is 4.6 MHz, yielding r eq = 129 /mi and A = 58 nm, for a result of K = 4.5 x 10~ 4 . The value 
of A relative to the crystal size is orders of magnitude smaller than in a quadrupole trap, meaning 
that interior ions in a multi-ion crystal should have correspondingly suppressed values of A. 

One remarkable feature of the data is the very high trap depth possible in a surface-electrode 
multipole trap. As an example, the depth of the trap used in Ref. ifPTI is approximately 170 meV, 
which is typical for millimeter- scale surface-electrode ion traps. For Trap A with r\ = 1 mm, 
we predict a depth of 50 eV when the trap is driven at 100 kHz with rf voltages of V\ = 100 V, 
V 2 = —13.66 V, and V 3 = 225.23 V. The reason why such a deep trap can be made is that the 
stability parameter 77 is not a constant within the trapping region, as is the case for quadrupole 
traps (in the region within which the electric potential is approximately quadratic), but depends in 
a nonlinear manner on r. 

Other example traps may be analyzed in a similar manner. Because of the somewhat large 
parameter space, we focus our attention. We choose a length scale of 1 mm and set the depth to 
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Trap 


Q/(2ir) [MHz] 


A [nm] 


V 


r eq [/^m] 




A 


4.6 


58 


1.0 x 10" 4 


129 


7.5 x 1(T 4 


B 


2.55 


121 


1.9 x 1(T 4 


194 


6.2 x 1(T 4 


C 


1.24 


919 


8.3 x 1(T 3 


213 


4.3 x 10~ 3 



TABLE IV: Table of values for the example traps. In all cases, the maximal rf voltage on any electrode is 
300 V, the length scale is set by r\ = 1 mm, and the trap depth is 100 meV. 



100 meV. We also select a maximum electrode voltage, choosing a common laboratory value of 



300 V. The values of r eq , D, and r] for one of two ions in each example trap are given in Table IV 
We do not further comment on the optimization of the ion height for a particular trap geometry; 
simple numerical experiments show that, at least for Trap A, z is already near the optimal value, 
given fixed values of the electrode radii. In this paper, we do not endeavor to present a truly 
"optimal" trap, since the criteria for what makes a trap optimal depend on the goals of a particular 
experiment. It may be possible, even just for the two-ion case, to improve on A for a fixed trap 
depth by searching the entire parameter space of r 2 , r 3 , and z . Once these are fixed, Vi and V3 
are uniquely determined, given an overall voltage scale V\ . 

We now consider higher-order traps. We observe that although the dependence of the pseu- 
dopotential on r and y is as expected, the trap depth falls precipitously as the order is increased. 
A plot illustrating this is given as Fig. |4} Strictly speaking, the traps should first be optimized 
in some way, but based on numerous numerical experiments (including different methods of op- 
timization), we think it likely that the unsealed depth is much smaller for all higher-order traps. 
The trap depth seen in Fig. [4] is not absolute; the real value depends on experimental parameters. 
Suppose we drive the trap at the highest voltage achievable, and then reduce £1 until the trap depth 
is some minimum acceptable value. A increases with just as D does. For Trap C, listed in Ta- 



ble IV A is roughly an order of magnitude higher than for a corresponding lower-order trap with 
the same trap depth and maximal voltage. K is also roughly an order of magnitude higher than the 
best lowest-order traps. Therefore we conjecture, that for all practical purposes, lowest-order mul- 
tipole traps are actually the best choice. Given sufficient optimization of the higher-order traps, 
however, this could turn out not to be true. 



Returning to the lowest-order trap, we again assume the trap potential is given by Eq. ( 10). To 



obtain the general formulas for 77 and A, we set Eq. ( [T2~| ) equal to the Coulomb force between the 
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FIG. 4: Logarithmic plots of ^(z) for Trap A for successive possible orders of p. The solid line represents 
the p = 6 trap, thick-dashed represents p = 10, and thin-dashed represents p = 14. The length scale for 
all traps is r% = 1 with zq = 1. The voltages for each trap are normalized such that the highest rf voltage 
applied to any trap electrode is 1 . Units are arbitrary. 

two ions. The equilibrium ion displacement r eq is then 



eq 



k e mVt 2 zl 



1/7 



(17) 



9Qc 2 V 2 _ 

where k e is the Coulomb constant: k e = 8.988 x 10 9 N-m 2 /C 2 . Taking the first and second 
derivatives of $(r) and using Eq. ^ and Eq. and substituting in our expression for r eq , we 
obtain the following expressions for r), A, and K: 



V = 2AQ 



V 3 k 2 



96 2 c r m 5 n 10 zl 2 



1/7 



(18) 



A = AQ 



c r Vk: 



96 3 m 4 fi 8 4j 



1/7 



(19) 



K = 4Q 



c 3 V 3 k 2 



96 2 m 5 Q 10 zl 2 



1/7 



(20) 



We can immediately check this model by examining the scaling of A with the frequency f2. 
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A fit to our data for Trap A shows that A oc f2~ L15 , and the predicted exponent in Eq. (19) is 
—8/7 ~ —1.14. We take this to mean that we can indeed approximate \I/(r) oc r 6 in the region 
of space occupied by two ions. Further confirmation has been obtained by comparing the values 



of A and r] computed by Eq. (19) and Eq. ( 18 ), respectively, to the values found by numerically 
computing E{r eq ) and VE(r eq ) according to the value of r eq found numerically. The values 
agree to within one percent, the error coming most likely from imperfect fits to the results of the 
numerical integration for \[>(r). 

Let us consider also the effect of the overall voltage scale: A oc V 1 ^ 7 /f2 8//7 . Using D oc V 2 /Vt 2 , 
we see that one can reduce the micromotion by raising both V and Vt so that their ratio remains 
constant. This is true in quadrupole traps as well, but with different powers in the scaling laws 
(rj oc V/VL 2 and A oc V 2 ^ 3 /f2 4 / 3 ). In practice, technical limitations will dominate the suppression 
of micromotion, even in the lowest-order traps. To give an example, suppose we want to use 
Trap A with z = 0.1 mm. With V x = 623 V and Q = 2vr x 286 MHz, A falls to 2.6 nm, with 
7] = 8.6 x 10~ 4 . These values, however, are already technically challenging. 

These scaling arguments highlight the marked differences between multipole and quadrupole 
traps. Whereas in a quadrupole trap the stability parameter q oc V/VL 2 and is a constant over the 
entire trap volume, the equivalent parameter r] depends strongly on the ion position, and in our 
examples here is far below 0.3. Furthermore, in a quadrupole trap, the pseudo-force T and the 
magnitude of the real electric field along r, E r , are both linear in the coordinate r. By contrast, T 
always scales up as a higher power of r than E, meaning that the ion is pinned to a significantly 
lower region of electric field than in a comparable quadrupole trap. This is, fundamentally, the 
reason for the desirable properties of ideal multipole traps: small values of 77, large possible values 
of D, and/ or the suppression of A. 



C. Multiple ions 

We now turn briefly to the question of the structure of multi-ion crystals in our multipole 
traps, and the related issues of stability and micromotion. We assume again that the temperature 
of the ions is zero, and find the positions that minimize the energy of the crystal. The results, 
not surprisingly, resemble the structure of crystals in linear multipole traps. In this article, we 
consider only the structures of small numbers of ions that form a ring in the plane parallel to 
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FIG. 5: Plots of crystal structures in Trap A with r\ = 1 mm, V\ = 100 V, and £l/(2ir) = 4.6 MHz, leading 
to a trap depth of 100 meV. Crystals of two, six, and twelve ions are represented, and the ion locations in 
each crystal are plotted with a "2," "6," and "12," respectively. In all cases, the displacement of each ion 
from z = zq is 17 fim or less, and thus the crystals are approximately planar. (This displacement increases 
with the number of ions.) 

the trap electrodes, assuming that ^(y) oc y 4 . We also specify the radial pseudopotential using 
Eq. (11). We present in Fig. [5]three representative crystal structures for Trap A. 



We may now easily calculate the micromotion amplitudes for ions in such 2-D crystals by 
calculating the magnitude of the electric field at the ion's location. We reference the crystals 
plotted in Fig. [5] For the two-ion crystal, A = 61 nm. For six ions, A = 142 nm, and for ions in 
the outer shell of the 12-ion crystal, A = 244 nm. However, the value of A for the interior ions 
is A = 13 nm. This raises the notion of using an outer shell of ions perhaps only for sympathetic 
cooling, while performing precision operations only on a ring of interior ions, if errors due to 
micromotion are a concern. 
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V. CONCLUSIONS 



We have shown, for the first time, a way to create an ideal multipole ion trap, in which ions 
are repelled from the rf null only by Coulomb repulsion and not by dc electric fields. Moreover, 
this is the first time that a surface-electrode multipole trap has been proposed. Although the basic 
procedure for solving for the voltages and/or electrode dimensions that null the quadratic terms in 
the pseudopotential is fairly straightforward, there is a large parameter space (including electrode 
widths and voltages, ion height, overall length and voltage scales, and drive frequency), which can 
make the route to optimization quite complicated. Here, we have tried to give a summary of some 
of the key features of these traps. 

Many questions remain. It would be advantageous to develop an automated routine that can 
optimize a trap under certain criteria, for certain numbers of ions. It is also not known whether it is 
possible to design higher-order SEMT's that are as deep as the lower-order versions while offering 
superior micromotion suppression. Additional concerns are the structure of multi-ion crystals in 
these traps, their motional frequencies and coupling rates, and their micromotion amplitudes. The 
computations of these quantities will require more sophisticated procedures and were beyond the 
scope of this paper. 

There are many potential applications of this work. The fact that no dc confinement is required 
means that micromotion, a source of systematic error (e.g. in atomic clocks [6|), can potentially 
be lower than in linear multipole traps. There is also the possibility of using SEMT's for the 
study of cold atom-ion collisions, since currently the temperature to which the mixture can be 
cooled is limited by the rf-induced heating of the ions 0. The utility of these traps for quantum 
computation and quantum simulation should also be considered. Since these traps are amenable 
to microfabrication, one could envision scaling down multipole ion traps, in a manner similar to 
that done with linear quadrupole ion traps in recent years. One implication is that ions can be 
placed near surfaces for quantum simulation using rf or microwave fields, in combination with 
magnetic field gradients |fT3ll25l . This is in addition to the benefit discussed in Sec. [TV] which is 
that micromotion amplitudes are suppressed, all else being equal, when the trap is scaled down. 
In the end, the applications seem numerous, and this list does not attempt to be exhaustive. 

We gratefully acknowledge funding from The Citadel and The Citadel Foundation and helpful 
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